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Summary of accomplishments 

During the course of the grant, we wrote and distributed about 12 reports and an 
equal number of journal papers supported fully or in part by this grant. The list of reports 
(title & abstract) and papers are given in Appendices A and B. 

This grant has indeed been instrumental in developing a robust hybrid Finite element 
method for the analysis of complex broadband antennas on doubly curved platforms. 
Previous to the grant, our capability was limited to simple printed patch antennas on mostly 
planar platforms. More specifically, 

• mixed element formulations were developed and new edge-based prisms were 
introduced. 

• these elements were important in permitting flexibility in geometry gridding for 
most antennas of interest. 

• new perfectly matched absorbers were introduced for mesh truncations associated 
with highly curved surfaces. 

• fast integral algorithms were introduced for boundary integral truncations reducing 
CPU time from 0(N 2 ) down to 0(N' 5 ) or less 

• frequency extrapolation schemes were developed for efficient broadband 
performance evaluations. This activity has been successfully continued by NASA 
researchers 

• computer codes were developed and extensively tested for several broadband 
configurations. These include FEMA-CYL, FEMA-PRISM and FEMA-TETRA 
written by L. Kempel, T. Ozdemir and J. Gong, respectively. 

• a new infinite balun feed was designed nearly constant impedance over the 800- 
3000 MHz operational band. 

• a complete slot spiral antenna was developed, fabricated and tested at NASA 
Langley. This new design is a culmination of the projects goals and integrates the 
computational and experimental efforts, this antenna design resulted in a U.S. 
patent and was revised three times to achieve the desired bandwidth and gain 
requirements from 800-3000 MHz. 


In the following we include; 

1. Appendix A: Reports Generated Under this Grant 

2. Appendix B: Refereed Journal Papers Supported Under this Grant 

3. Appmdix C: Abstracts of Reports Generated Under this Grant 

4. Summary viewgraphs on the computational capabilities of the developed hybrid 
FEM incorporating fast algorithms 

5. An update on the latest analysis, measurement and performance of the patented slot 
spiral 

6. A review of the hybrid finite element method, much of which was developed under 
this grant 
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Forward 


This grant began in January 1993 and ended December 31, 1997. The goals of the 
grant were to develop analysis and design techniques for conformal broadband slot or 
printed antennas mounted on doubly curved surfaces. 

The uniqueness of the project can be summarized as follows: 

• a first implementation of conformal antennas on doubly curved platforms 
(cylindrical, spherical or any other non-planar configuration). 

• development of suitable finite element formulation which combined various 
elements and domain decompositions for treating the fine details of printed 
broadband antennas such as spirals and log periodics. 

• investigation of efficient mesh truncation schemes for enclosing the finite element 
domain. New absorbers and boundary integral enclosures with fast integral 
algorithms were examined 

• the ultimate design, fabrication and testing of a broadband conformal antenna which 
covers all mobile communication needs from 800 MHz up to 3 GHz. 
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APPENDIX A 

Reports Generated Under This Grant 


Report # 

Title 

Authors 

030601- I-T 
(July 1993) 

Radiation and scattering by cavity-backed 
antennas on a circular cylinder 

Leo C. Kempel, 
John L. Volakis 

03060 1-2-T 
(September 1993) 

Scattering by a groove in an impedance plane 

Sunil Bindiganavale, 
John L. Volakis 

03060 1-3-T 
(January 1994) 

Scattering by cavity-backed antennas on a 
circular cylinder 

Leo C. Kempel, 
John L. Volakis 

03060 1-4-T 
(January 1995) 

Simulation of thin slot spirals and dual circular 
patch antennas using the finite element method 
with mixed elements 

J. Gong, J. L. Volakis, 
M.W. Numberger 

03060 1-5-T 
(July 1995) 

Triangular prisms for edge-based vector finite 
element analysis of conformal antennas 

J. Gong, T. Ozdemir, 
J. L. Volakis 

03060 1-6-T 
(August 1995) 

An new planar feed for slot spiral antennas 

M.W. Nurnberger, 
J. L. Volakis 

03060 1-7-T 
(January 1996) 

Efficient finite element simulation of slot 
antennas using prismatic elements 

J. Gong, 

J. L. Volakis 

030601 -8-T 
(July 1996) 

An AWE implementation for electromagnetic 
FEM analysis 

J. Gong, 

J. L. Volakis 

03060 1-9-T 
(August 1996) 

A planar slot spiral for conformal vehicle 
applications 

M.W. Numberger, 

J. L. Volakis, D. T. Fralick, 
F. B. Beck 

030601- 10-T 
(October 1996) 

Robust hybrid finite element methods for 
antennas and microwave circuits 

J. Gong, 

J. L. Volakis 

030601-1 1-T 
(August 1997) 

Fast memory-saving hybrid algorithms for 
dectromagnetic scattering and radiation 

S. S. Bindiganavale, 
J. L. Volakis 

030601- 12-T t 

(Janaury 1998) t 

\WE method in frequency domain 
dectromagnetics 

E. Erdemli, C. Reddy, 
J. Volakis 


3 


















APPENDIX B 

Refereed Journal Papers Supported Under this Grant 

1. L C. Kempel and J.L. Volakis, "Scattering by cavity-backed antennas on a circular 
cylinder," IEEE Trcms. Antennas Propagat. , Vol. 42 . No. 9, pp. 1268-1279,1994. 

2. J.L. Volakis and L. C. Kempel, "Electromagnetics. Computational methods and 
considerations," IEEE Computational Science and Engineering Mag., Vol, 2 . No. 1, 
pp. 42-57, Spring 1995. 


3. L.C. Kempel, J.L. Volakis and R. Sliva, "Radiation by cavity-backed antennas on a 
circular cylinder," IEE Proceedings-H, pp 233-239, 1995. 

4 J.L. Volakis and A. Chatterjee, " A selective review of the finite element-ABC and 
the finite element-boundary integral methods for electromagnetic scattering," (invited 
paper) Annales de Telecom. , Vol, 50 . No. 5-6, pp. 499-509, May- June 1995. 

5. J. Gong and J.L. Volakis, “An efficient and accurate model of the coax cable feeding 
structure for FEM simulations,” IEEE Trans. Antennas and Propagat Vol 43 nn 
1474-1478, Dec. 1995. 

^ J Gong, J.L. Volakis and H T.G. Wang, “Efficient finite element simulation of slot 
antennas using prismatic elements,” Radio Sci. , Vol. 3 1, No 6 pp 1837-1844 Nov- 
Dec. 1996 

7. J.Gong and J.L. Volakis, “An AWE Implementation for Electromagnetic FEM 
Analysis,” IEE Electronics Letters, Vol. 32, No.24, pp.2216-2217, 21st November 
1996. 

8. M. Numberger and J.L. Volakis, “A new planar feed for slot spiral antennas,” I EEE 
Trcms. Antennas and Propagat., Vol. 44, pp. 130-131, Jan. 1996. 

9. S. Bindiganavale and J. L. Volakis, "Scattering by narrow gaps in an impedance 
plan e," Radio Sci.., pp 401-408, March-April 1996. 


10. J.L. Volakis, J. Gong and T. Ozdemir, “Applications of the FEM to Conformal 
Antennas, Ch. 13 in the book Finite Element Software for Microwave Engineering 
ed. by Itoh, Silvester and Pelosi (Wiley, 1996), pp. 3 13-345 


11. T. Ozdemir, M.W. Numberger, J.L. Volakis, R. Kipp and J. Berne, “A hybridization 
of finite element and high frequency methods for pattern prediction of antennas on 
aircraft structures,” IEEE Antennas and Propagat. Soc. Mag, Vol. 38 No 3 pp 28- 
38, June 1996. 


12 J.L Volakis, T. Ozdemir and J. Gong “Hybrid Finite Element methodologies for 
antennas and scattering,” IEEE Trans. Antenna Propagat., pp 493-507, March 1997 

13 T. Ozdemir and J. L. Volakis, “Triangular Prisms for Edge-Based vector finite 
elements analysis ,” IEEE Trans. Antennas Propagat., pp 788-797, May 1997 

^ ^ Li, P. Papalambros, J.L. Volakis, ‘Designing broadband patch antennas using the 
sequential quadratic programming method,” I EKE Trans. Antennas and Propaeat 
pp. 1689-1692, Nov. 1997. ’* 

15. S. Bindiganavali, J. L. Volakis and H. Anastassiu, “Scattering from plates containing 
small features using the adaptive integral method(AIM),” submitted to IEEE Trans. 
Antennas and Propagat. 

16. S. Bindiganavale and J.L. Volakis, "Radiation and Scattering from Cavity-Backed Slot 
Antennas Using the Finite Element-Adaptive Integral Method," submitted to IEEE 
Trans. Antenna Propagat 

AMTA Papers 

L M. Numberger, J. Volakis, D. Fralick and F.B. Beck, ” A planar slot spiral for 
conformal vehicle applications,” 1997 American Measurement techniques 
Assoc. (AMTA) Symposium, Seattle, Digest pp. 121-126. 


2. M.W. Numberger and J.L. Volakis, "A planar slot spiral for mobile communications," 
19th American Measurement Techniques Assoc. (AMTA), Symposium Digest pp 
Boston, MA, 1998. 


SEVERAL papers were also presented at the following Meetings: 

1 Antennas and Propagation Symposia (held annually in North America) 

2. Applied Computational Electromagnetics Society 

3. International URSI Meetings 
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APPENDIX C 


Abstracts of Reports Generated Under 
this Grant 



030601- 1 -T 


Radiation and Scattering by Cavity- backed 
Antennas on a Circular Cylinder 

Leo C. Kemp el and John L. Volakis 
July 21, 1993 


Abstract 

Conformal arrays are popular antennas for aircraft and missile 
platforms due to their inherent low weight and drag properties. How- 
ever, to date there has been a dearth of rigorous analytical or nu- 
merical solutions to aid the designer. In fact, it has been common 
practice to use limited measurements and planar approximations in 
designing such non-planar antennas. In this report, we extend the fi- 
nite element-boundary integral method to scattering and radiation by 
cavity-backed structures in an infinite, metallic cylinder. In particular, 
we discuss the formulation specifics such as weight functions, dyadic 
Green’s function, implementation details and particular difficulties in- 
herent to cylindrical structures. Special care is taken to ensure that 
the resulting computer program has low memory demand and mini- 
mal computational requirements. In this report, both scattering and 
radiation parameters are computed and validated as much as possible. 
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03060 1 -2-T 


Scattering by a Groove in an Impedance 

Plane 

Sunil S.Bindiganavale and John L.Volakis 

Radiation Laboratory 
Department of Electrical Engineering 
and Computer Science 
The University of Michigan 
Ann Arbor, Michigan 48109-2122. 

September 24, 1993 


Abstract 

An analysis of two-dimensional scattering from a narrow groove 
in an impedance plane is presented. The groove is represented by 
a impedance surface and hence the problem reduces to that of scat- 
tering from an impedance strip in an otherwise uniform impedance 
plane. On the basis of this model, appropriate integral equations are 
constructed using a form of the impedance plane Green’s functions in- 
volving rapidly convergent integrals. The integral equations are solved 
by introducing a single basis representation of the equivalent current 
on the narrow impedance insert. Both transverse electric (TE) and 
transverse magnetic (TM) polarizations are treated. The resulting 
solution is validated by comparison with results from the standard 
boundary integral method (BIM) and a high frequency solution. It is 
found that the presented solution for narrow impedance inserts can 
be used in conjunction with the high frequency solution for the char- 
acterization of impedance inserts of any given width. 
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03060 13- T 


Scattering by Cavity-backed Antennas on a 
Circular Cylinder 

Leo C. Kempcl and John L. Volakis 
Radiation Laboratory 
1301 Beal Ave. 

Ann Arbor. MI 48109-2122 

November 15. 1993 


Abstract 

Conformal arrays are popular antennas for aircraft, spacecraft and 
land vehicle platforms due to their inherent low weight and drag prop- 
erties. However, to date there lias been a dearth of rigorous analyt- 
ical or numerical solutions to aid the designer. In fact, it has been 
common practice to use limited measurements and planar approxi- 
mations in designing such non -planar antennas. In this paper, wo 
extend the finite element-boundary integral method to scattering by 
cavity-backed structures in an infinite, metallic cylinder. In particular, 
we discuss the formulation specifics such as weight functions, dyadic 
Green's function, implementation details and particular difficulties in- 
herent to cylindrical structures. Special care is taken to ensure that 
the resulting computer program has low memory demand and mini- 
mal computational requirements. Scattering results are presented and 
validated as much as possible. 
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03060 1-4 1 ' 


Preface 

This semi-annual report describes our progress up to mid-January 1995. The report contains five 
sections all dealing with the modeling of spiral and patch antennas recessed in metallic platforms. 
Of significance is the development of decomposition schemes which separate the different regions 
of the antenna volume. These schemes are used in modeling the spiral and the dual circular patch 
antennas. In the case of the spiral, a mixed (edge and node) element formulation is developed to 
avoid use of small elements in modeling the slot. This also allows for use of the simpler wire 
green s function in modeling the radiated fields. Overall, the resulting modeling approach leads to 
substantial simplifications and is perhaps the most efficient within the context of hybrid FEM 
formulations. The implementation of this formulation is nearly complete but we already fabricated 
and tested an archimedean slot spiral operating at a center frequency of 1.1 GHz (see Section 2 of 
the report). The main goal of the experiment was to examine the feasibility of a new feed design 
for the slot spiral antenna. This feed design consisted of a microstrip line below the metallic 
section of the slot spiral but, unfortunately, our tirst design trial required improvements to bring 
the axial ratio close to unity. A new design is currently under fabrication and w ith the help of the 
new software, it is our hope that the feed design issue will be resolved over the next two months. 

Substantial effort was devoted over the past six months in improving the feed model in the 
context of the FEM. This is addressed in Sections 3 and 4 w'here we describe simple, but very 
successful implementations of the coax probe and slot-coupled microstrip line feed. In the latter 
case, the feed network is modeled using brick elements whereas the antenna volume is modeled 
using tetrahedrals. The two computational regions were then coupled by enforcing potential 
continuity (a new approach) across the slot joining the two regions. A design of a dual stacked 
circular patch configuration is included which demonstrates the validity of the potential continuity 
condition. Also, as noted in Section 3, the proposed coax cable model yields input impedance 
results which overlay the measured data. 

Finally, Section 5 of the report presents an innovative scheme for truncating finite element 
meshes. This is an alternative to the absorbing boundary conditions and employs a perfectly 
matched anisotropic absorber. It was already tested for terminating transmission lines in 
microwave circuits and we are in the process of implementing it in our FEM antenna codes. 

In the next few months our goal is to complete the spiral antenna code and employ it to examine 
several spiral antenna configurations operating at different bands, including that from 944- 
1256 MHz as discussed in the proposal. The code will be also enhanced with several features 
relating to feed modeling and mesh truncation. Finally, additional fabrications and antenna 
measurements will be made 



030601-5 T 


Triangular prisms for od go- based vector finite 
element analysis of conformal antennas 

1 Ozdcmir and .]. L. \ olakis 
Radiation Laboratory 

Department of Electrical Engineering and Computer Science 

University of Michigan 
Ann Arbor, Michigan 48109-2122 

June 5. 1995 

Abstract 

I his paper deals with I he derivation and validation of edge- based shape fnn< 
lions for the distorted I riangular prism. Although I he let rahedron is often the 
element of choice lor volume tessellation, mesh generation using tetrahedra i>. 
cumbersome and ( 1*1 intensive. On the other hand, the distorted triangnlai 
prism allows for meshes which arc unstructured in two dimensions and sirin' 
lined in the third dimension. This leads to substantial simplifications in tin 
meslnng algorithm and many printed antenna and microwave circuit groin.- 
I tics can be easily tessellated using such a mesh. The new edge-based shape 
functions presented in tins paper are validated by eigenvalue, radiation and 
scattering, and input impedance computations involving conformal antennas 
on planar and cylindrical platforms. Also presented are results showing t he 
effect of platform curvature on resonance behavior of such antennas. 
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A New Planar Feed for 
Slot Spiral Antennas 


M. \V. Xuriiborgcr and J. L. Volakis 

Radiation Laboratory 
Depart incut of Electrical Engineering 


and Computer Science 
rhe University of Michigan 


Ann Arbor, MI 48109 2122 


Abstract 

I his report presents a new planar, wideband feed network for a 
slot spiral antenna, and the subsequent design and performance of a 
\ lit antenna utilizing this feed design. Both input impedance and 
radiation pattern measurements are presented to demonstrate the per- 
formance and usefulness of this feed. Almost all previous designs have 
utilized wire spirals, requiring bulky, non-planar feeds with separate 
bain ns, and large absorbing cavities. The presented slot spiral an- 
tenna feed integrates the balim into the structure of the slut spiral 
antenna, making the antenna and feed planar. This greath simplifies 
the design and construction of the antenna, in addition to providing 
repeatable accuracy ft also allows the use of a very shallow reflecting 
cavity for conformal applications Finally, this feeding approach now 
makes man\ of the known miniaturization techniques viable options. 
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03060 1- 7- T 


Efficient Finite Element. Simulation of 
Slot Antennas Using Prismatic Elements 


Jian Gong and John L. Volakis 
radiation laboratory 
University of Michigan 
Ann Arbor, Ml 48109-2122 
jxgum@umich.edu volakis@umich.edu 


Abstract 

A hybrid finite element - boundary integral (FE-BI) simulation 
technique is discussed to treat narrow slot antennas etched on a planar 
platform. Specifically, the prismatic elements are used to reduce the 
redundant sampling rates and ease the mesh generation process. Nu- 
merical results for an antenna slot and frequency selective surfaces are 
presented to demonstrate the validity and capability of the technique. 


1 Introduction 

It has been reported that a hybrid finite element boundary integral tech- 
nique [Jin, et. al. 1991, Silvester and Pelosi (199-1)] can be employed for 
characterizing conformal antennas of arbitrary shape [Gong, et. al. (1994)]. 
Indeed, planar/non-planar, rectangular/non-rectangular designs, ring slot 
or spiral sloi antennas with probe, coax cable or mlcrostrip line feeds can be 
simulated with this formulation. This is because of the geometrical adapt- 
ability of tetrahedral elements used for the implementation. However, in 
practice, certain configurations require extremely high sampling rates due 
to the presence of fine geometrical details. Among them are a variety of 
slot antennas (spirals, rings, slot spirals, cross slots, log-periodic slots, etc.), 
where the slot width is much smaller than the other dimensions (cavity di- 
ameter or inter-distance of slots). In these cases, the mesh is extremely 
dense (with over 50, 100 or even higher samples per wavelength), whereas 
typical discretizations involve only 10 20 elements per wavelength. This 

1 


13 
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An AWE Implementation for 
Electromagnetic FEM Analysis 

•liaii Gong and John L. Volakis 
Radiation Laboratory 
University of Michigan 
Ann Arbor. MI 481098-2122 
.jxguinfUiniich.odu volakis'Ojuinich.odn 


Abstract 

A It hough full wave elec t romagnetic systems are large and mm 
bersonie to solve, typically only a few parameters, such as input 
impedance. S parameters, and far field pattern, are needed by the* 
designer or analyst. A reduced order modeling of these parameters is 
therefore an important consideration in minimizing the CPE require 
menls. I lie Asymptotic Waveform Evaluation! AWE) method is one 
approach u> construct a reduced order model of the input impedance 
or other useful electromagnetic parameters. We demonstrate its appli- 
cation and validity when used in conjunction with the finite element 
method to simulate full wave electromagnetic problems. 


1 Introduction 

The method ol Asymptotic Waveform Evaluation (AWE) provides a rerlucoci- 
order model ol a linear system and has already been successfully used in VLSI 
and circuit analysis to approximate the transfer function associated with a 
given set ol pot t s/variables in circuit networks [1, 2, 3], The basic idea of 
the method is to develop an approximate transfer function of a given linear 
system i from a limited set of spectral solutions. Typically, a Pade expansion 
of the transfer lime (ion is postulated whose coefficients are then determined 
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A PLANAR SLOT SPIRAL FOR 
CONFORMAL VEHICLE APPLICATIONS 

M.W. Nurnberger*, J.L. Volakis’, D.T. Fralick', F.B. Beck 1 


* Radiation Laboratory 
Dept, of Elect. Engin. & Comp. Sci. 
University of Michigan, Ann Arbor, MI 48109-2122 

f NASA-Langley Research Center 
Hampton, VA 23681 


ABSTRACT 


A slot spiral antenna and its associated feed are presented for conformal mounting on a variety of 
an , air, and sea vehicles. The inherent broadband behavior and good pattern coverage of the spiral 
antenna is exploited for the integration of multiple frequencies, and thus multiple transmitting and 
receiving apertures, into one compact, planar antenna. The feasibility of the broadband slot spiral 
antenna relies on the use of an equally broadband, balanced, planar, and non-intrusive feed struc- 
ture The design of the slot spiral, its feeding structure, and the reflecting cavity are discussed with 
emphasis on the experimental validation and construction of the antenna. 

Keywords: Spiral Antennas, Planar, Conformal, Automotive.-Antenna Measurements 


1.0 INTRODUCTION 

Spiral antennas are particularly known for their ability to produce very wide-band, almost perfect 
circularly-polarized radiation over their full coverage region. Because of this polarization diveristy 
and broad spatial and frequency coverage, many applications exist, ranging from military surveilancc 
bCM, and ECCM to numerous commercial and private uses, including the consolidation of multiple’ 
low gain communications antennas on moving vehicles. 

For the typical wire spiral antenna, the performance advantages mentioned above come at the 
pnee of size and complexity. While the radiating elements of a wire spiral may be planar, the feed 
network and balun structures generally are not, and combine to add weight, depth, and significant 
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030601- I 0-T 


Robust Hybrid Finite Element Methods 
for Antennas and Microwave Circuits 


Introduction 


One of the primary goals in this dissertation is concerned with the development 
of robust hybrid finite element-boundary integral (FE-BI) techniques for modeling 
and design of conformal antennas of arbitrary shape. Both the finite element and 
integral equation methods will be first overviewed in this chapter with an emphasis 
on recently developed hybrid FE-BI methodologies for antennas, microwave and 
millimeter wave applications. The structure of the dissertation is then outlined. We 
conclude the chapter with discussions of certain fundamental concepts and methods 
in electromagnetics, which are important to this study. 
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ABSTRACT 

FAST MEMORY-SAVING HYBRID ALGORITHMS FOR 
ELECTROMAGNETIC SCATTERING AND RADIATION 

The main theme in this dissertation is examination of techniques for alleviating the 
excessive memory requirement and execution time associated with implementations 
of popular frequency domain numerical techniques for electromagnetic scattering 
and radiation. Of particular interest is the reduction of computational complexity 
associated with integral equation techniques such as the moment method. 

In this dissertation, we develop an algorithm for directly generating sparse ma- 
trices from the original fully populated moment method interaction matrix and sub- 
sequently refining the solution iteratively. However, this technique is still O(N^) 
and yields more accurate results for relatively smooth scatterers. We were the first 
to perform an error analysis of the Fast Multipole Method (FMM) and based on 
given error criteria, we set guidelines for choosing the various parameters affecting 
the speed and solution accuracy. Scattering from a realistic size aircraft nose radome 
and the interactions between a radome and the antenna at its base was computed 
with a version of the FMM. For the first time, fast integral methods are considered for 
efficient solution of hybrid finite element boundary integral methods, thus enhancing 
the attractiveness of these techniques to model large, composite bodies. 

The Adaptive Integral Method (AIM) is ideal for modeling bodies with intricate 
construction and high discretization rates since it translates the original problem 
grid to a new coarser regular grid, enabling use of Toeplitz matrix techniques for fast 
matrix-vector product computation. We demonstrate special features that make this 
technique eminently suitable for planar scatterers and cavity backed-antennas. 
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AWE Technique in Frequency Domain Electromagnetics 


Yunus E. Erdemli’, C.J.Reddy* and John L. Volakis’ 


Abstract: 

This paper presents fast radar cross section (RCS) computations using the 
Asymptotic Waveform Evaluation (AWE) technique in conjunction with Method of 
Moments (MoM) and hybrid Finite Element (FEM/MoM) implementations. In its 
traditional form, AWE constructs a reduced-order model of a given linear system by 
Taylor series expansion with respect to specific values of the system parameters 
(frequency, angle, etc.). Thus, AWE permits the prediction of the frequency response 
from a few frequency calculations. In this paper we modify AWE to instead employ a 
rational function (Pade approximation) representation for the system parameters. Using a 
Pade rational function instead of a Taylor series, the accuracy of the analysis is increased 
to a wider frequency range. AWE is also extended to allow monostatic RCS pattern 
prediction using again a few pattern values, thus eliminating a need to resolve the system 
when an iterative solver is employed. To demonstrate these extensions of AWE, 
numerical examples of three-dimensional metallic bodies and cavity-backed apertures are 
considered. 


Y. E. Erdemli and J.L. Volakis are with Radiation Laboratory, University of Michigan, Ann Arbor MI, USA. 
* C. J. Reddy is with Department of Electrical Engineering, Hampton University, Hampton VA 23668, USA. 
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Hybrid Finite Element Methods for 
Conformal Antenna Simulations 
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Finite Element Method-History 

FEM is a mature method 
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Many FEM success for Antennas and Circuits 
Accurate Feed Modeling 
Can Avoid Green’s Functions 
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•Number ot Unknowns is reduced by a factor of 5 or more 
•For Narrow slots BI matrix is small 

•FFT can be used to speed-up computation of matrix vector products 
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TRIANGULAR CAVITY SLOT 
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Other Antennas Analyzed With 
Hybrid Finite Element Techniques 
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Slot Spiral Performance and Design for 
Personal Communications 
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using periodic array 
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memory by a factor of 9 
CPU by a factor of 3 
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0(n 2 ) vs. 0(n* 5 ) Memory Locations 
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Design Optimization Approach 
in Conjunction with FEM 
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Slot Fed Dual Patch Antenna: Design for a Specified 
Bandwidth and Frequ.Using SQP: Varying 3 Parameters 
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A PLANAR SLOT SPIRAL FOR 
MULTI FUNCTION COMMUNICATIONS APERTURES 

M.W. Numberger and J.L. Volakis 


Radiation Laboratory 
Dept, of Elect. Engin. & Comp. Sci. 
University of Michigan, Ann Arbor, MI 48109-2122 


ABSTRACT 

A slot spiral antenna and its associated feed are presented 
for conformal mounting on a variety of land, air, and sea 
vehicles. By exploiting the inherent broadband behavior, 
good pattern coverage and polarization diversity of the spi- 
ral antenna, a conformal antenna which can be concur- 
rently used for cellular, digital personal communications 
(PCS), global positioning (GPS) and intelligent vehicle 
highway systems (WHS) as well as wireless LAN net- 
works has been developed. A key requirement for achiev- 
ing such broadband behavior (800-3000MHz) is the avail- 
ability of a broadband planar feed and balun. Such a feed 
was proposed last year by the authors. However, addi- 
tional design improvements were found to be necessary to 
achieve satisfactory pattern and gain performance. Among 
them were a broadband termination for the spiral arms and 
the suppression of cavity and waveguide modes. Both of 
these improvements played a critical role in achieving ac- 
ceptable performance over the 800-3000 MHz bandwidth. 
After a general description of the slot spiral antenna and 
the above modifications, this paper presents a comparison 
of the performance before and after the modifications. 

Keywords: Spiral Antennas, Planar, Conformal, Automo- 
tive, Antenna Measurements, Broadband, Circular Polar- 
ization 

1.0 INTRODUCTION 

Spiral antennas are particularly known for their ability 
to produce very wide-band, almost perfect circularly- 
polarized radiation over their full coverage region. Be- 
cause of this polarization diversity and their inherent broad 


spatial and frequency coverage, a wide range of applica- 
tions exist. These range from military surveillance, ECM, 
and ECCM to numerous commercial and private uses, in- 
cluding the consolidation of multiple low gain communi- 
cations antennas on moving vehicles. 

For the typical wire or microstrip spiral antenna, the 
performance advantages mentioned above come at the 
price of size and complexity. While the radiating elements 
of a wire spiral may be planar, the feed network and balun 
structures generally are not, and combine to add weight, 
depth, and significant complexity to the system. Further- 
more, because the spiral antenna radiates bi-directionally, 
a broadband absorbing cavity is generally used to elimi- 
nate the radiation in one direction, adding even more depth 
to the antenna. While some designs exist [1] that integrate 
the feed and balun into the cavity and reduce the com- 
plexity somewhat, the absorbing cavity is still at least a 
quarter- wavelength deep at the lowest frequency of opera- 
tion, adding significant thickness to the antenna [2]. 

A slot spiral is not burdened with most of these diffi- 
culties. As was demonstrated previously [3], the balun and 
feed structure may be integrated into the planar radiating 
structure. This greatly simplifies the construction and in- 
creases the accuracy of the antenna by allowing standard 
printed circuit techniques to be used throughout the en- 
tire fabrication process. Also, while the slot spiral also 
radiates bi-directionally, a deep, absorbing cavity is not 
necessary to force uni-directional radiation. In fact, for 
conformal mounting, a very shallow reflecting cavity is 
sufficient, and also serves to increase the gain somewhat. 

Of particular interest in this paper is the design of a 
slot spiral antenna which occupies minimum real estate 
(small thickness and diameter) and operates at least up to 
3 GHz. More specifically, we are interested in integrat- 
ing the transmit/receive functions of the following systems 
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into a single aperture: 

• Analog and Digital Cellular Telephone (825-950 
1700- 1900 MHz) 

• Global Positioning System (GPS) (1227.6 and 
1575.42 MHz) 

• Personal Communications System (PCS) (1850- 
2000 MHz) 

• Intelligent Vehicle Highway System (IVHS) 

• IEEE 802. 1 1 Wireless LAN (2400-2500 MHz) 

Hence, a minimum bandwidth of 3: 1 is required to suc- 
cessfully integrate these functions. Although the earlier 
design could easily handle this bandwidth, to maintain a 
small antenna diameter (less than 6 inches), we found it 
necessary to develop both better spiral arm terminations 
and new cavity mode suppression techniques. In the fol- 
lowing, we first present a description of the overall antenna 
and then proceed to discuss its performance after the afore- 
mentioned improvements. 


2.0 SLOT SPIRAL DESIGN 

The design procedure for the slot spiral antenna assem- 
bly is somewhat long and involved, and requires the care- 
ful choice of many inter-related design parameters. In 
this section we summarize the concerns in the design of 
the slot spiral, and discuss the more important design 
parameters. The design parameters are: spiral growth 
rate, a; number of turns, N; microstrip line characteristic 
impedance(s), Zm. ; slot line characteristic impedance(s), 
Z s .; substrate characteristics: e r , tan <5, and metaliza- 
tion technique; thickness of the substrate, f; length and 
impedance of the tuning stub, l stub & Z Stub \ cavity depth, 
h ; cavity wall material; slot line termination material and 
technique; and the specific geometry of the microstrip-slot 
line transition. 

As with most antennas, the design procedure of the 
slot spiral is an iterative one. However, because the in- 
dividual parameters are highly dependent on each other, 
the design procedure is particularly involved. It is perhaps 
easiest initially to choose the growth rate, and then to care- 
fully inspect the initial design for inconsistencies. For the 
following discussion, the reader is referred to Figures 1 - 
4 for clarification. 

Physically, the crucial section of the antenna is a sin- 
gle printed circular substrate whose top and bottom views 
are displayed in Figures 1 and 2. On one side of the 20 


mil substrate (see Figure I) is etched the slot spiral radi- 
ating element. The other side (see Figure 2) of the sub- 
strate incorporates the infinite balun microstrip feed. This 
conductor operates as a microstrip transmission line using 
the other side as a ground plane, i.e. the planar conductor 
forming the slot spiral. As seen, the microstrip line spi- 
rals without interfering with the slot radiation to the center 
of the antenna. At the center, the microstrip is allowed to 
cross the slot spiral and is either extended approximately 
a quarter wavelength beyond the feed point or is brought 
through the substrate with a via, where it is shorted to 
the opposite side of the slot. When no ground plane is 
employed, the antenna radiates bi-directionally, with sym- 
metric patterns on both sides. A ground plane is necessary 
in practice for conformal mounting on vehicles (airframe 
skin, car roof, boat skin etc.) to serve both as a shield and, 
more importantly, to eliminate the need for absorber back- 
ing. The cavity can also be used to improve the gain by as 
much as 3 dB. 

Because of the slot operation, the slot spiral radiation 
can be interpreted in terms of equivalent magnetic cur- 
rents whose radiation is enhanced with the presence of a 
ground plane in the immediate proximity of the slot. Thus, 
a ground plane or cavity (see Figure 3) in conjunction with 
a slot spiral allows for very thin antenna designs. In addi- 
tion, it provides shielding and nearly doubles the gain over 
a very wide bandwidth (larger than that of the antenna). 
We should remark that the side or vertical walls of such 
cavities are typically terminated with some absorbing ma- 
terial to suppress cavity resonances, but as found in our in- 
vestigation, additional cavity treatment may be necessary 
to realize the broadband antenna characteristics. More- 
over, as noted above, the slot arms should be terminated 
properly to suppress the traveling wave reflections which 
interfere both with the polarization purity of the antenna 
and its bandwidth. Without any loading, it would be nec- 
essary to resort to a much larger antenna size which is, 
of course, highly undesirable for the intended commercial 
uses of the antenna. 

Having discussed the overall antenna geometry and its 
crucial parameters, we now proceed with the specifics of 
the design. In the past [3] we presented a design of the 
slot spiral and the infinite balun microstrip feed. However, 
it was later found that additional modifications were nec- 
essary, including a new slot arm termination and cavity 
and waveguide mode suppression. Simulation and mea- 
sured results displaying the characteristics of the antenna 
after the implementation of the tapered resistive slot ter- 
mination are shown later. Also, we discuss the importance 
of the cavity design for proper antenna operation over the 
entire bandwidth. 

Let us first review the basic slot spiral and feed de- 


52 


sign. In the 500 MHz to 4 GHz region, we found that a 
growth rate of 1 10.9 mils/rad to be a good initial choice. 
The initial characteristic impedance of the microslrip line 
was chosen to be 50ft. As can be realized from the above 
antenna description, the actual feed is contrived through 
the use of a microstrip to slot line transition. Here, the mi- 
crostrip line views the slotline as a pair of shunt branches. 
Thus, for the transition to operate most effectively, the 
characteristic impedance of the microstrip line is chosen 
to be half that of the slot line, yielding a perfect match at 
the feed. Based on the dielectric constant and thickness of 
the substrate, and on fabrication concerns, a certain range 
of slot line impedances will be available. For a dielec- 
tric constant of 3.38 and a thickness of 20 mils, a 1 10ft 
slot line, which has a width of approximately 35 mils, was 
used. This required a microstrip characteristic impedance 
of 55ft at the transition, and yielded a microstrip ground 
plane width of approximately 313 mils. 

As the microstrip feed spirals to the slot spiral cen- 
ter, of concern is the coupling between the microstrip feed 
line and the slotline, the latter being the primary radia- 
tor. Any coupling from the microstrip feed to the slotline 
will have adverse effects on the pattern shape and polar- 
ization purity of the antenna. To minimize coupling, it 
is best to maximize the ratio of microstrip line to ground 
plane width. Many design issues can be traded off in this 
case - the most important being the growth rate, substrate 
thickness and dielectric constant, slotline impedance and 
the microstrip line characteristics. In our design approach, 
we first chose to modify just the microstrip line width, and 
adjust the others if necessary later. To increase the width 
ratio mentioned above, the impedance of the intermediate 
microstrip line was increased until the width ratio was ap- 
proximately doubled. We also employed Klopfenstein ta- 
pers [4] in both cases to transition from the 50ft outer feed 
impedance to the 70ft intermediate microstrip impedance, 
and then back to the 55ft microstrip line impedance at the 
center of the spiral. In this manner, the impedance changes 
do not hinder the broadband behavior of the antenna and 
feed, thus maintaining the slot spiral's broadband behav- 
ior. 

3.0 SLOT SPIRAL RESISTIVE 
TERMINATION 

Several measurements of the original 6 inch diameter and 
0.2 inch thick (with ground plane) spiral antenna were col- 
lected at the NASA Langley facilities. The antenna was 
tested in isolation, and a 1 /4-scale version was also tested 
on the roof of a 1 /4-scale minivan, with the goal of assess- 
ing its bandwidth, pattern coverage, CP purity, gain and ef- 
ficiency. Specifically, measurements were performed from 


500 MHz to 5 GHz, with a representative pattern given in 
Figure 5. We show the pattern of the antenna in the mid- 
dle of the roof of the minivan at the full scale frequency of 
1350 MHz, obtained with a spinning linear source. Thus, 
the rapid oscillations observed in the pattern are a measure 
of the axial ratio of the antenna at this frequency. In this 
case, zero degrees refers to radiation normal to the spiral, 
and 90 degrees corresponds to radiation in the plane of 
the antenna. As expected, the antenna gain is maximum 
along the broadside and drops somewhat in the plane of 
the antenna. At 70 degrees, the gain is down approxi- 
mately 5 dB, and at 80 degrees, about 10 dB. The axial 
ratio is also expected to worsen (higher oscillations are as- 
sociated with worse axial ratio) at large angles from nor- 
mal but for this frequency, the axial ratio needs improve- 
ment even at broadside (zero degrees). This is true for the 
patterns of the isolated antenna (not shown) and for the 
antenna mounted in the middle of the minivan's roof. Ac- 
tually, towards the lower frequencies, the antenna's overall 
performance deteriorates substantially. This deterioration 
was traced to traveling wave reflections from the slot spi- 
ral arm termination. To maintain the same antennU size 
while minimizing these reflections and improving its low- 
frequency performance, it was necessary to apply some 
type of tapered absorber termination over the ends of the 
spiral. 

Initially, tapered applications of foam and ferrite- 
loaded rubber absorbers were used, but these techniques, 
although they demonstrated the importance of the termina- 
tion, did not lead to a repeatable design for manufacturing 
purposes. Instead, we developed a technique using resis- 
tive paint applied in a tapered fashion over the slots. The 
compactness of the resistive paint was found to be very at- 
tractive, and the overall shape and thickness of the termi- 
nation was also consistently reproducible in manufactur- 
ing. To verify the design of the resistive paint termination, 
the commercial simulation package IE3D was employed, 
permitting the modeling of the substrate, cavity and the 
entire spiral configuration, including the resistive paint in 
the slots, via the moment method. A display of the ge- 
ometry used by IE3D is shown in Figure 4 and includes 
the microstrip feed, the conductors forming the slot spi- 
ral, and the resistive termination of the slots. As is the 
case with most moment method codes, the conductors are 
broken into triangular elements and a current expansion 
is introduced over all elements. The latter are found by 
enforcing the integral equation with the multi-layer sub- 
strate Green's function on all conductor segments to form 
a linear system of equations. 

Samples of the simulation results for the antenna using 
the resistive taper are shown in Figures 6 and 7, and corre- 
spond to frequencies of 800 MHz and 3 GHz. Calculations 
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were also performed at many intermediate frequencies, but 
aren’t shown here. In Figures 6 and 7 we display the left 
and right circularly polarized patterns for two orthogonal 
cuts of the antenna (<£ = 0 and $ = 90). As usual, the 
elevation angle is 8, and 8 — 0 refers to broadside radia- 
tion. Three patterns are included in each figure: the first 
refers to operation of the isolated antenna with no cavity 
and without the resistive termination; the second pattern 
refers to the antenna after the resistive taper has been ap- 
plied, but still with no cavity. The third pattern is for the 
final antenna configuration, which includes both the resis- 
tive termination as well as the cavity. 

It is clear from the patterns in Figures 6 and 7 that at 
cellular frequencies, the antenna is nearly useless without 
the use of the resistive termination (bad axial ratio, non- 
uniform radiation, etc.). However, when the termination 
is applied, the expected performance of the slot spiral is 
recovered with a gain of nearly 5 dB tc at broadside and 
nearly 20 dB of polarization isolation. The calculated an- 
tenna efficiency in this case is slightly larger than 70%. 
When the cavity is introduced, the pattern shape in general 
is maintained, but the effects on the polarization purity are 
dependent on frequency, as might be expected. The effi- 
ciency also drops substantially. These difficulties are due 
to cavity and waveguide mode excitations which can in 
part be suppressed by introducing absorber at the cavity 
walls. However, these absorbers create other problems of 
their own, not to mention a certain loss of efficiency. Al- 
ternatively, other approaches can be employed to suppress 
the cavity modes while minimizing losses in efficiency. 
Because of space limitations, these remedies will be dis- 
cussed at the conference. As seen from the patterns in 
Figure 7, at higher frequencies, while the axial ratio per- 
formance is quite good, pattern symmetry and coverage 
still are negatively effected by the cavity, as is the effi- 
ciency. 


4.0 CONCLUSIONS 

In this paper we presented a design for a slot spiral an- 
tenna which has the bandwidth and pattern coverage nec- 
essary to integrate cellular, GPS, PCS, and wireless lan 
system functionality into a single aperture. The key fea- 
tures of this antenna are its small thickness (less than 0.2 
inches) and small diameter (less than 6 inches), thus per- 
mitting conformal and inconspicuous mounting on mobile 


platforms. At the same time it integrates several communi- 
cation and navigation systems into a single platform, thus 
reducing the overall cost of the integrated system while 
maintaining versatility. 

The proposed antenna achieves its small thickness by 
using an innovative design which incorporates a planar in- 
finite balun whose bandwidth is maintained over the en- 
tire antenna operation. It also incorporates novel tapered 
resistive arm terminations and cavity mode suppression 
schemes to allow the realization of the full potential of the 
antenna. The effectiveness of the resistive taper and the 
overall operation of the antenna was verified with a CAD 
simulator based on the method of moments, and measure- 
ments were collected at specific frequencies. Some mea- 
surements were also collected with a 1 /4-scale antenna 
mounted on the roof of a 1/4-scale minivan. In all cases, 
the resistive arm termination proved essential to the proper 
operation of the antenna, particularly at lower frequencies. 
Cavity mode suppression was also viewed as essential. To- 
gether, these remedies should permit the effective use of 
slot spiral antennas which, in spite of their attractiveness, 
have so far been of limited usefulness. 
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Figure 1 : Underside of Slot Spiral Antenna, showing 
slot radiator. 


Figure 2: Top view of Slot Spiral Antenna, showing in- 
finite balun microstrip feed. 



Figure 3: Bottom view of Slot Spiral Antenna, showing 
absorber walls and reflecting cavity bottom. 



Figure 4: Geometrical model of antenna for numerical 
analysis. 
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Figure 5: Spinning linear plot of spiral antenna on the roof of a minivan, showing poor axial ratio. 
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Introduction 

Wide band properties and circular polarization characteristics of spiral an- 
tennas make them quite attractive for mobile and modern telecommunication 
systems. However, their deployment on conformal metallic platforms has been 
impeded by adverse effects from the ground plane or the cavity. Many studies 
have been carried out to investigate ways of successfully incorporating ground 
planes into the spiral design. Nakano et. al. [1] and Wang, et. al. [2], to name 
a few, showed that metal backing is compatible with the spiral modes through 
their studies of mostly wire and printed spirals. 

In this study, we show that cavity-backed slot spirals (rather than printed) 
are particularly attractive for conformal applications. Since slot radiation is 
enhanced by the presence of a shallow metal-backed cavity, very thin (less than 
a ^ thick) antennas can be constructed without compromising the spiral’s 
broadband performance. 

A specific design proposed in [3] is considered below for analysis and im- 
provement. 


Design 

Figure 1 shows the geometry of the cavity-backed spiral. The spiral is desired 
to have a gain of at least 3dB with less than 2dB axial ratio from 800 to 
3000 MHz [3]. Slot radiating elements have been chosen because they tolerate 
very thin reflecting cavities, which produce a unidirectional pattern and a 3dB 
increase in gain. The cavity depth was set at 0.635cm thus varying from A/59 
to A/16 over the operational band. This thickness is, of course, much thinner 
than typical absorbing cavities. The antenna is fed with a microstrip feed 
with an integrated balun [4], which is planar and easy to build. The slot 
termination is important in achieving circular polarization (low axial ratio) 
and a resistive thin-film termination was developed to allow arbitrary taper 
functions for optimum matching. The size of the antenna determines the low 
frequency performance and the number of turns is chosen such that the outer 
circumference is larger than the largest operational wavelength. 
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Finite element analysis with mixed-elements 

To analyze the complex spiral geometry and its cavity, the field is formulated 
via the Finite Element Method (FEM) in conjunction with edge-based ele- 
ments [5]. The Boundarj' Integral (BI) technique has been used to terminate 
the mesh at the cavity aperture. The implementation for the BI follow closely 
that in [6]. The uniqueness of our analysis approach lies more with the meshing 
than the implementation of the field solution. A mix of hexahedra and prisms 
[5] were used to model the volume resulting in a combination of quadrilaterals 
and triangles on the aperture. This permitted an accurate and physical mod- 
eling of the aperture (which is perforated with the thin slot spiral) without 
resorting to oversampling. The particular mesh shown in Figure 2(a) results in 
a total of 1307 unknowns, a rather small size. One drawback is the elongated 
elements in the slot which worsen the condition number of the system matrix. 
Nevertheless, no difficulty was encountered in achieving tolerances of less than 
one percent using a biconjugate gradient solver. 

In the simulations given next, the antenna was excited by placing a current 
source across the slot at the center of the antenna. The input impedance was 
then obtained by dividing the computed voltage (across the slot) at the feed 
location by the amplitude of the excitation current. Also, the resistive film 
taper was modeled with equally spaced equivalent lumped resistors whose 
values were chosen to dissipate an equal amount of energy. 


Antenna performance 

First, we comment on the fundamental operations of the slot spiral. Fig- 
ure 2(b) shows the computed variation of the voltage (magnetic current) along 
the spiral arms. The real and imaginary parts are 90 degrees out of phase indi- 
cating a traveling wave (also implying that the resistive termination eliminates 
reflections from the spiral ends and this is necessary for low axial ratio). Fur- 
ther, the separation of the wave crests is one free-space wavelength indicating 
a pure TEM wave. Basically the fields tightly couple to the narrow slot and 
the propagation constant is not altered by the presence of the ground plane. 
However, the ground plane is responsible for the voltage lobes toward the spiral 
ends, a feature also present with printed spirals [1], Note that, at 1900MHz, 
the voltage distribution is approximately the same as that of a center-fed 1A 
linear dipole. 

Next, let us look at the broadband behavior of the antenna. The results 
in Figure 3 show that the bandwidth increases as more turns are added or 
the expansion ratio is increased. The low frequency gain improves while the 
axial ratio remains at acceptable levels. Results also show that increasing the 
expansion ratio a is slightly better than adding more turns (same expansion of 
the antenna size) in terms of gain performance. The decrease in gain toward 
the low end of the frequency band is due to the resistive termination of the 
spiral. As the frequency drops, the radiating slots move toward the rim of the 
spiral and thus more of the energy is dissipated in resistive termination. Since 
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Figure 1: (a) Cavity-backed Archimedean slot spiral. 

the resistive termination is essential for low axial ratio, this calls for a compro- 
mise between the gain and the pattern. Figure 3(b) shows that the axial ratio 
can be maintained less than 2dB using the resistive termination but achieving 
gain above zero dB requires increase in diameter to about 20cm. Note that 
in all simulations, the cavity thickness was not altered and the microwave ab- 
sorbing ring was placed at a radial distance of 1.2cm from the slot radiator. 
Bringing it closer increased the overall gain but deteriorated the axial ratio. 
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Figure 2: (a) Surface mesh comprised of quadrilateral and triangular patches, 
(b) Voltage along the spiral arms for 800 and 1900 MHz. 
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Figure 3: Slot spiral performance as a function of expansion ratio a and number 
of turns: (a) Gain, (b) axial ratio. 
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Abstract 

This paper is an overview of the finite element method (FEM) as applied to 
electromagnetic scattering and radiation problems. A brief review of the method- 
ol °g.V is given with particular emphasis on new developments over the past five 
years relating to feed modeling, parallelization and mesh truncation. New applica- 
tions which illustrate the method’s capabilities, versatility and utility for general 
purpose application are discussed. 


1 Introduction 

Over the past 10 years we have witnessed an increasing reliance on computational meth- 
ods for the characterization of electromagnetic problems. Although traditional integral 
equation methods continue to be used for many applications, one can safely state that 
in recent years the greatest progress in computational electromagnetics has been in the 
development and application of partial differential equation(PDE) methods such as the 
finite difference-time domain and finite element (FEM) methods, including hybridiza- 
tions of these with integral equation and high frequency techniques. The major reasons 
for the increasing reliance on PDE methods stem from their inherent geometrical adapt- 
ability, low O(N) memory demand and their capability to model heterogeneous (isotropic 
or anisotropic) geometries. These attributes are essential in developing general-purpose 
analysis/design codes for electromagnetic scattering, antennas, microwave circuits and 
biomedical applications. For example, modern aircraft geometries contain large non- 
metallic or composite material sections and antenna geometries may involve many lay- 
ers of materials, including complex microwave circuit feeding networks. Although, the 
moment method continues to remain the most accurate and efficient approach for sub- 
wavelength size bodies of simple geometries, PDE methods and hybrid versions of these 
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have shown a much greater promise for large scale simulations without placing restric- 
tions on the geometry and material composition of the structure. 

This paper is a selected review of hybrid finite element methods and their applica- 
tions to scattering and antenna analysis. We begin by first introducing the mathematical 
basics of the method without reference to a specific applications and in a manner that 
identifies the method's inherent advantages in handling arbitrary geometrical configura- 
tions which may incorporate impedance boundary conditions and anisotropic materials. 
In this section we also identify the key challenges associated with the implementation of 
the method such as mesh truncation and feed modeling for antenna applications. Section 
3 of the paper reviews the various mesh truncation schemes to be employed on the outer 
mesh surface for the unique solution of the vector wave equation. Absorbing boundary 
conditions, integral equations and artificial absorbers are discussed, all leading to differ- 
ent versions of hybrid FEM methodologies, and we comment on their accuracy and ease 
of implementation. Section 4 presents several approaches for antenna feed modeling in 
the context of the FEM, including coaxial as well as aperture coupled feeds. The next 
section is devoted to parallelization issues specific to finite element codes. We give per- 
formances of typical FEM codes and provide storage and implementation guidelines for 
maximizing code performance on parallel computing platforms. The final section of the 

paper gives some representative applications of the method to scattering and antenna 
problems. 


2 Theory 

2.1 FEM Formulation 

Consider the antenna and scattering configurations shown in Figure 1. In the case of 
a scatterer, the entire computational domain is enclosed by a fictitious surface S Q that 
may encompass a variety of composite/dielectric volumes as well as metallic, impedance 
and resistive surfaces. The antenna geometry is assumed to be recessed in some doubly 
curved surface. In this case, the bounding surface S„ may either be located as shown or 
can be coincident with the antenna aperture. As usual, the recessed cavity is intended to 
house the radiating elements and their feeding structure such as coaxial cables, striplines, 
microstrip lines or aperture coupled feeds. The cavity may encompass any inhomoge- 
neous or anisotropic material including resistive cards, lumped or distributed loads and 
so on. 

The goal with any finite element formulation is to obtain the solution of the vector 
wave equation 

V x p r • (V x E)] - k% • E = - jKZJi - V x (p; 1 • M,-) (1) 

in which E represents the total electric field, (c r ,/7 r ) denote the relative permittivity and 
permeability tensors of the computational domain, ko is the free space wavenumber and 
Z 0 is the free space intrinsic impedance. In addition, J, and M, denote the impressed 
electric and magnetic sources, respectively, and represent the excitation for the antenna 
problem. As is well known, the standard finite element (FE) solution scheme is to 
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Figure 1: Computational domains for FEM analysis, (a) the various regions enclosed by 
S o, (b) typical tetrahedral mesh, (c) computational regions for antenna analysis 


consider the weak form of (1), instead of solving it directly. This can be achieved by 
extremizing the pertinent functional [1] 


F(E) = 


+ 

+ 

+ 


5 JJJ {(V X E) • py • (V X E)j - e 0 (l ■ E) • e} dv 
JJJ ® ' [j^oZoJ t + V x (jI T 1 * M,j] dv 

jk Q Z 0 If E • (H x n ) ds 
j J Sq+$J 


-jkoZti 


{Us, A 


x E) ■ (h x E) ds + 


d 

Zls 


IIl EEdv } 


( 2 ) 


where R denotes the resistivity or impedance of the surface S r , S s stands for the junction 
opening to possible guided feeding structures and H is the corresponding magnetic field 
on S 0 and Sj whose outer normal is given by h. Also, V s is the source volume and Vi is 
the volume occupied by the load Zl, whose length and cross section are given by d and 
s, respectively. It is noted that S r must be closed when it satisfies an opaque impedance 
boundary condition but can be open (i.e. a finite plate) if it satisfies the penetrable 
resistive sheet condition [2] 


fix fix E = -Rn x (H + - H“) 

where denote the fields on the two sides of the surface S r . As seen from (2), in spite 
of the discontinuity in the magnetic field, no special care is required for the evaluation of 
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tlir integral over S T - The explicit knowledge of H is, however, required over the surface 
.So and Sf (referred to as mesh truncation surfaces) for the* unique solution of E. 

As an alternative to (2). we can instead begin with (1 ) by weighting it with a testing 
function T. Subsequently, application of the divergence theorem yields the residual 


< R, T > = 


+ 

+ 

+ 


/ flv { (V X T) ‘ W ' (V X E) ] “ k ° l * • E • T } dv 

J j j v T • {jk o Z 0 J, + V x (^T 1 • M,-)} dV 
jk o Z 0 ff T • (H x h)dS 


j h Q Z 0 


Us, ^ 


X T ) • (n x E )dS + 


-irsUL, 


T • EdV 


( 3 ) 


This is typically referred to as the weak form of the wave equation, whereas (2) is referred 
to as the variational form. It will be seen later that when set equal to zero, (3) yields 
the same system of equations as those deduced from the functional (2). Therefore both 
methods are equivalent. 

When dealing with large computational domains as is often the case in scattering 
applications, it is instructive to work with the scattered rather than the total field, and 
also distinguish the air and dielectric regions. This approach results in reduced errors 
within the computational domain and to proceed we introduce the definitions 


E = E !nc + E sca< , H = H"* c + H ica ‘ (4) 

where E tnc denotes the excitation field impinging from the exterior of S 0 and E 3ca * is the 
scattered field. With these substitutions, the functional F takes the form 


F(E scai ) = 


+ 

+ 

+ 

+ 

+ 

+ 


^ IIIv X E *“ f ) ' ( V x ESCOi ) ~ ' E scat } dv 

J fj v E scai • [jkoZoJ, + V x (f; 1 • M,-)] dV 
jk 0 Z 0 J j 


Sq + Sj 


(H x n) dS 


E scai . E sca( dV 


l -jk 0 Z 0 | JJs i(n x E scat ) • (n x E scat ) dS + ^ J jj 
jk 0 Z 0 [fj s i(n x E scat ) • (n x E inc ) dS+J-jjJ^ E scat ■ E inc dV J 

\ JJ J {v x E scat • (^T 1 • V x E' nc ) - klf r ■ E' ne • E scat 

m 


V x E 


tnc I — 


(/^r 


. \7 x E scai 


'} dV (5) 

) - jfcg? r ■ E scai • E' nc ) dV + /(E inc ) 


in which / (E mc ) is a collection of terms involving only the incident field E ,nc . Its specific 
form is of no consequence since it will be eliminated in the subsequent minimization of 
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F. Most importantly, we observe in (0) the presence of integrals over the volume \ ' d 
occupied by the dielectrics and over the surface S r . The introduction of these integrals is 
required when working with the scattered field because E iVf1 ' is discontinuous at bound- 
aries separating media with different permeabilities [4], Also note that the presence of 
the last two integrals is necessitated because E" ,c satisfies the wave equation only in free 
space, and that these two volume integrals become equal when the dielectric anisotropy 
reduces to isotropy. 

The functionals (2) and (6) or the weighted residual (3) can be cast into a discrete 
system of equations for the solution of E. To accomplish this, it is first necessary to 
subdivide the volume as a collection of small elements such as those shown in Figure 2 
[3]. Within each volume element, the field can then be expanded as 
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Cylindrical Shell 


Figure 2: Illustration of the various elements for tessellating three'dimensional volumes 


M =#of edges 

E'= £ EfWl = {W'\ T {E‘} (6) 

k=l 

in which W£ are referred to as the edge-based shape or basis functions of the eth element 
in the computational domain. In contrast to the traditional node-based shape functions, 
the edge-based shape functions are better suited [5] for simulating three dimensional elec- 
tromagnetic fields at corners and edges. Moreover, edge-based shape functions overcome 
difficulties associated with spurious resonances [6]. They were proposed by Whitney [7] 
over 35 years ago and were revived in the 1980s [8], [9], Their representation is different 
for each element but have the common properties [10] of being divergence free (i.e., 

5 


65 



V ' Wl ; = 0 within ,llf ‘ ) and normalized in such a manner so that the expansion 

coefficients E\ : represent the average field value across the At h edge of the all element. 
One disadvantage of the edge-based elements is that they increase the unknown count. 
However, this is balanced by the increased sparsity of the resulting stiffness matrix. A 
detailed mathematical presentation of the edge-based shape functions for various two 
and three dimensional elements (bricks, hexahedra and tetrahedra) [1 1]-[14] is given by 
Graglia ft. al. [15] in this issue of the Transactions. Linear as well as higher order 

(curvilinear) elements are discussed [16]-[19]. 


2.2 Discretization and System Assembly 

Once the computational domain has been tessellated using appropriate elements and 
shape functions, the discretization of (2), (3) or (6) proceeds by introducing the expan- 
sion E or E 5Ca< = Y, e =! E e , with E e as given by (6) and N e denoting the number of 
elements comprising the computational domain. For the functional (2), the system of 
equations is constructed by setting (Rayleigh- Ritz procedure) 


dF( E) 

m 


= 0 , 


e — 1,2, A r e ; k = 1,2 , M 


( 7 ) 


whereas in the case of (3) Galerkin’s approach is used by setting T = W£. Regardless 
of the procedure, the resulting system will be of the form 


+ £[*•]{£*} + £>•'} + ^{CU = 0 (S) 

e=1 *=1 e=l p -1 

in which the brackets denote square matrices and the curly braces refer to columns 
Among the various new parameters, [A e ] is referred to as the element matrix and results 
from the discretization of the first volume integral in (2), (3) or (6); N s is equal to the 
aggregate of the surface elements (quadrilaterals for hexahedra or triangles for tetrahedra 
and prisms) used for the tessellation of S 0 , Sj and S r \ the column {A' e } results from the 
discretization of the source integral over V,] and [tf s ] is the matrix associated with the 
third and fourth surface integrals in (2), (3) or (6). Finally, N p denotes the aggregate 
of the surface and volume elements over S T , Ve and V d and the column {C p } results 
from the discretization of the integrals in (6) involving the external incident field Eb 
Basically, { C p } provides the external excitations (scattering problem) whereas {A' e } is 

the corresponding excitation column for internal sources as is the case with antenna 
problems. 

The entries of the element matrix [A e ] are specific to the choice of the shape functions 
and are compactly given by 

M*i = Jfl {[DW'm\-'{DW'f - k][W'm\W'] T } dV (9) 
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with 


[D ir] r 


£ - jg: {"7} r 

i-AK) 7 -t,{^} t 


( 10 ) 


and l t denotes the volume of the discrete element while the subscripts (r.j/.c) in (10) 
imply the (x.y, z) components of the shape function W f . If lumped loads are present 
(i.e., in the presence of a volume integral over If), the diagonal entries of [,4 e ] are simply 
modified with the addition of the term jk 0 Z 0 d 2 /Z L . Perfectly conducting posts located 
at the k th edge are handled by setting the corresponding total fields equal to zero and 
rearranging the final matrix as discussed later. 

The internal source excitation column is simply given by 



and the entries of the corresponding excitation column due to the incident field with the 
presence of dielectric materials are 



in which W p implies the representation of the shape function over the required p ih 
surface or volume element. 

The specification of the matrix [fi s ] can not be completed without first introducing 
some a priori relationship between the E and H fields on S 0 . This relationship (or 
boundary condition on S 0 ) is referred to as the mesh termination condition and its form 
depends on the physics of the problem. For example, in the case of problems where the 
entire computational domain is enclosed by a perfect electric or magnetic conductor, the 
unknown column {£ s } is eliminated from {E}. For scattering and radiation problems, 
the mesh termination condition should be such that the artificial boundary S 0 appears 
transparent to all waves incident from the interior of the computational domain. Clearly, 
if S 0 is placed far enough from the scatterer or radiator, the simple Sommerfeld radiation 
condition provides the appropriate relationship between E and H. However, this is not 
practical since it will yield large computational domains. To bring S 0 as close as possible 
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to the scattering or radiating surface, more sophisticated boundary conditions must be 
introduced on .s' 0 . These mesh termination conditions are critical to the accuracy and 
efficiency of t he formula! ion and are some of t he major hot t lenecks in t he implement at ion 
of the FEM. Regardless of the employed mesh termination approach, after carrying out 

the sums in (8), a process referred to as matrix assembly, the resulting matrix system 
will be of the general form 



In this, [E l } denotes the field unknowns within the volume enclosed by S 0 + S/ whereas 
{E } represents the corresponding unknowns on the boundaries S 0 and/or Sj. The ex- 
citation column {b v } results from the assembly of {A' e } and similarly {6 5 } is associated 
with {CP}. The matrix [A] is very sparse (9 to 30 non-zero entries per row) and this 
is a major characteristic and advantage of FEM. By using special storage schemes and 
solvers suitable for sparse systems, the CPU and memory requirements are maintained 
at O(N) and scalability can be attained on multiprocessor platforms. 

The boundary matrix [Q] in (12) is associated with the boundary fields {E B } and its 
specific form is determined by the employed mesh termination condition on S 0 . Over the 
past five years, much research on FEM has concentrated on the development of mesh 
truncation schemes which minimize the computational burden without compromising 
accuracy. As will be seen later, this compromise is difficult to attain and is subject 
to the computational problem being considered. For the purpose of our discussion, we 
will subdivide the various mesh termination schemes into three categories: (1) exact 
boundary conditions, (2) absorbing boundary conditions (or ABCs), and (3) artificial 
absorbers. Exact boundary conditions provide an integral relation between the electric 
and magnetic fields, and because they are exact, they permit the placement of S 0 very 
near or exactly on the surface of the scatterer or radiator. The resulting formulation is 
referred to as the finite element-boundary integral (FE-BI) method and combines the 
adaptability of the FEM and the rigor of the boundary integral methods. However, it 
yields a fully populated matrix [Q] and as a result, the FE-BI method is associated with 
higher computational demands even though only a small part of the overall system is full 
(see Figure 3a). To alleviate the higher computational demands of the rigorous FE-BI 
method, ABCs or artificial absorbers (AAs) can instead be used to terminate the mesh. 
Both of these are approximate mesh termination methods, but lead to completely sparse 
and scalable systems of equations. Basically, the sub-matrix [Q] is eliminated with the 
application of ABCs or AAs (see Figure 3b). In the case of ABCs, a local boundary con- 
dition in the form of a differential equation is applied on S 0 to relate E and H so that S 0 
appears as transparent as possible to the incident fields from the interior. The resulting 
method will be referred to as the finite element-ABC (FE-ABC) method and has so far 
been the primary method for general-purpose scattering computations. Finally, the use 
of artificial absorbers (including perfectly matched layers or PMLs) [20] have recently 
gained major attention because of their potential for greater accuracy and inherent im- 
plementation simplicity. In the context of the finite element-artificial absorber (FE-AA) 
method, the mesh is terminated by using a material absorber (typically non-realizable 
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in practice) to absorb the outgoing waves and suppress backward reflections. Below, we 
briefly discuss the specifics for each of the three mesh termination schemes. 


3 Mesh Termination Schemes 

3.1 Finite Element-Boundary Integral Method 

The FE-BI method was introduced in the early seventies [21], [22] as a natural exten- 
sion of the FEM for modeling unbounded problems. However, because of its larger 
computational requirements, the method did not enjoy a widespread application to elec- 
tromagnetics until the late eighties [23, 24], In accordance with the FE-BI method, the 

relation between E and H on S 0 is provided by the Stratton-Chu integral equation (or 
its dual) 

h x H(r) = ri x H inc (r) 

+ fix II {[n' x H(r)] x V'G 0 (r, r') + jk 0 Y 0 [n' x E(r)]G 0 (r,r') (13) 

J J S Q 

+ 2 ^v , -[n' x E(r')]VG 0 (r,r')} dS' 

where r and r denote the observation and integration points, respectively, and G 0 ( r, F) = 
exp(-jk 0 \r - r'|/(47r|r - r'|) is the free-space Green’s function. The above is the’ most 
general form of the boundary integral and places no restrictions on the shape of S a with 
the exception that it must be closed as shown in Figure 1(a). Provided S 0 is placed 
just above the outer boundary of the scatterer or radiator, any material composition 
which is interior to S 0 can be handled with ease using the FEM without relation to the 
boundary integral. This form of the FE-BI was used by Yuan [25] and later by Angelini 
tt. al. [26], Antilla and Alexopoulos [27], and Soudais [28] to model three dimensional 
scatterers with anisotropic treatments. The method has also been successfully used for 
biomedical simulations [29], [30]. 

The discretization of (14) follows by introducing the expansion (over the s th surface 
element) 


of edges 

E(r)= £ £JS:(R)-[5*] t {£‘} (14) 

fc=l 

with a similar expansion for the H field (if necessary). The surface shape functions 
S s (r) can be set equal to W e (r) when the latter are evaluated on the surface of the 
volume element. For the tetrahedron, n x S k are simply equal to the traditional basis 
functions given by Rao et. al. [31]. However, as pointed out in [32], S s can be chosen 
independently of W e provided care is exercised when (14) is substituted/combined with 
(2)-(6). Regardless, when (14) is substituted into (3) or (6) after discretization we will 
get the partly sparse, partly full system [33] given in (12) and illustrated in Figure 3(a). 
Finally, we point out that since S 0 is a closed surface, the final system is not immune to 
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false interior resonances clue to the boundary integral. In this case, the combined field 
formulation [34] or the simpler complexification scheme [35] may be employed. 

When S 0 coincides with an aperture in a ground plane (see Figure 1(b)). the inte- 
gral equation (14) can be simplified substantially [36]. Specifically the integral relation 
between E and H on the aperture reduces to 


nxH = ft x H 3 ° -f 


JL* x 


G(r, r') [E(r) x ft] dS 


(15) 


where G is the dyadic Green’s function of the second kind, with ft x V x G = 0 on S a 
and for planar platforms it reduces to , 


G = ~j~ ( I + ^VV 


z. 


kl 


e -jfco|r-r'| 

47r jr — r'| 


(16) 


In this, I is the unit dyad and the factor of 2 is due to image theory. Also, is 
equal to the sum of the incident and reflected fields (in the absence of the aperture) for 
scattering or zero for antenna parameter computations. For non-planar S Q , H 3 ° is equal 
to the field scattered by an otherwise smooth surface again with the aperture removed. 
In this case, the Green s function should also be modified accordingly with respect to 
the non-planar platform. One of the first implementations of the FE-BI for radiation 
and scattering from rectangular apertures/antenna recessed in a ground plane was given 
by Jin and Volakis [14], [33], [36] and was later extended to antenna analysis on planar 
[37], [38], [39] and cylindrical [40] platforms. 

The FE-BI method is particularly attractive in terms of CPU and memory require- 
ments when the [Q] matrix is Toeplitz and can therefore be cast in circulant form [41], [42], 
[43], [44], In this case, the entire system can be solved using an iterative solver [45] in con- 
junction with the FFT to reduce the CPU requirements down to 0(N se log N se ) for the 
boundary integral sub-matrix. The FFT is simply used to carry out the matrix-vector 
product [Q] { E B } within the iterative solver. For example, if the symmetric biconjugate 
gradient (BCG) method [45] is employed to solve (12) and rectangular elements are used, 
the storage requirement is only 4A r ee + 8 N ae , where N ee and N se denote the number of 
edges within the computational volume and on S 0 , respectively. For triangular surface 
elements, the storage requirement is about 4.5 times larger due to the presence of the 
diagonal edges across the quadrilaterals. Whether the full matrix \Q] is Toeplitz or not 
depends on the shape of the surface S 0 and on the uniformity of the mesh. It can be 
shown that [Q] will be Toeplitz for planar and cylindrical surfaces provided the surface 
grid is uniform. The Toeplitz form of [£?] has been exploited with great success and sys- 
tems with more than 0.5 million unknowns spanning apertures of 15A x 15A have been 
solved accurately on a desktop workstation. Convergence of the iterative solver is also 
quite impressive for these systems and as demonstrated by Ozdemir and Volakis [46] (see 
Figure 4) the FE-BI method is the method of choice for these applications. For example, 
we observed that a 180,000 unknowns system converged in 57 iterations and for another 
system of 25,000 unknowns convergence was achieved after 66 iterations with an average 
CPU time of 2 sec/iteration on an HP9000/750 rated at 24.7 MFlops. Because of these 
impressive performances, it is advantageous to transform [£] to a Toeplitz matrix even 
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ulirn the surface grid is not uniform. To do so. a uniform grid can he superimposed onto 
the non-uniform mesh (see Figure 5) and the edge fields on the 1\vo grids can he related 
via a connectivity matrix [39]. In this manner, non-rectangular antenna elements and 
apertures can he modeled with the full geometrical adaptability of the finite element 
method and without compromise in accuracy. It is important to note that similar trans- 
formation matrices can be exploited for decomposing computational domains as done, 
for example, in [47]. 



Number of Array Elements in x direction 

Figure 4: Comparison of the convergence behavior/CPU requirements among the FE-BI, 
FE-ABC and FE-AA methods 

When S 0 is not planar, the boundary integral matrix will inevitably cause large CPU 
and storage requirements for large N ae to the point where FE-ABC or FE-AA methods 
become more attractive at the expense of accuracy. Recently, though, techniques have 
been proposed which show promise in reducing the computational requirements of the 
matrix-vector product [Q] {E 5 }. Among them, the fast multipole method [48], [49] can 
reduce the operation count in carrying out the matrix- vector product from O ( N] e ) down 
to O {N} e 5 ) and reductions down to O (A^ 33 ) have also been proposed. The main idea of 
the FMM [49] is to subdivide the surface S 0 into groups, each containing approximately 
M se ~ \fN7 e unknowns. By rewriting the free space Green’s function as an expansion [50] 
or by introducing its far field approximation [51], it can be shown that the interactions 
between the unknowns within the groups can be carried out in O (N se M se ) operations 
whereas the interactions between groups can be carried out in 0(N]JM se ) operations. 
The sum of these two operation counts yields a total operation count of 0(N}*) for the 
boundary integral matrix when M se fa s/NZ, and this must be added to the O (N ee ) 
operation count associated with the sparse matrix [.4]. Details for the implementation 
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circular p,n<4. 



Figure 5: Overlay of a uniform grid on an unstructured mesh for implementation of the 
FFT to perform the matrix-vector products. 


of the FMM are given in [48], [49] and when combined with the FEM, we can refer to 
it as the FE-FMM method [52]. A comparison of the FE-BI (with LU and conjugate 
gradient-FFT solvers) and the FE-FMM methods for the calculation of the scattering 
by a large groove in a ground plane is tabulated in Table 1. The groove was 50A long 
and 0.35A deep (all metallic) and filled with a dielectric having e r = 4 and fj, r = 1. By 
using a sampling rate of about 15 edges per wavelength, the system unknowns were 
Nee = 5256 and N se = 750. From the data in Table 1, it is clear that the FE-FMM is 
more than 3 times faster and requires 3-5 times less memory than the FE-BI when LU 
decomposition is used for solving the FE-BI system. Because of the grouping/averaging 
of the surface elements, the FE-FMM exhibits certain errors which are functions of many 
parameters [53] and these errors must be kept in mind as the system size is increased. It 
is expected that the FMM will exhibit greater speed-ups as the system size is increased. 
Nevertheless, the special implementation of the FE-BI using the FFT is still by far the 
most accurate approach. 

We close this section by noting that another approach for treating the boundary 
integral matrix- vector products is via the adaptive integral method (AIM) [54], This 
approach relies on the introduction of multipole expansions (similar to Taylor series ex- 
pansions) to replace the sources on S 0 by equivalent point sources placed on rectangular 
grids. In this manner the CPU requirements can be reduced by making use of three di- 
mensional FFT routines. It can again be shown that AIM reduces the boundary matrix 
CPU requirements down to O (N} e 5 ) or even down to 0 (A^ 33 ). 
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Groove Width 

Total Unknown*. 

BI Unknowns 

CPI TiTnt Jor HI i .\finittr 5. seconds ! 

Fh-Bl (CO) 

i-T-FMM 

FF.-BI (C'GFFT) 

25 A 

2f>31 

375 

{SA 8) 

(3.26) 

(1.11) 

35 A 

3681 

525 

(16.31) 

r (5.55) 

(3,2-1) 

50 A 

5256 1 

750 

(45.1) 

(14.31) [ 

(■5.40) 


Groove Width 

Storage of BI (KB) 

FE-BI (CG) 

FE-FMM 

to 

Ol 

1072 

277 

35A 

2102 

458 

50A 

4291 

784 



RMS error (dB) 

Groove Width 

FE-FMM 

25 A 

1.12 

35A I 

1.2 

50A 

1.36 


Table 1: CPU Times, storage requirement and error 

3.2 Finite Element-Absorbing Boundary Condition Method 

The goal with any ABC is to eliminate backward reflections from S Q , and a variety of 
these boundary conditions have been proposed over the years beginning with those of 
Bayliss and Turkel [55], and Engquist and Majda [56]. More recently, other ABCs such 
as those by Webb and Kanellopoulos [57], and Chatterjee and Volakis [4] (see also Senior 
and Volakis [2]) have been proposed. All ABCs provide an approximate relation between 
the E and H fields on the surface So which in most cases is derived by assuming a field 
expansion in inverse powers of r, where r is the radial distance from the center of S 0 . If 
the ABC annihilates the first (2m + 1) inverse powers of r, it is then referred to as an 
mth order ABC. The second order ABC derived by Chatterjee and Volakis [4] is given 
by 


jk 0 Z 0 (H scat x n) = or • E* cai + /3 • [V x n(V x E scaJ ) n ] + 7 • V*(V • E° cat ), ( 17 ) 
where 

Q — + Ot 2 i 2 i 2 

P = + t 2 t 2 ) 

T = -i U^o + 3k to — — ^ 2«i)f 1 t 1 + -7£-(jk 0 + 3«: m - — 2K 2 )i 2 t 2 

K m jk Q Km ' 

Q 1 — K g T D(jko — /Cl) T + AC m A/c] 

«2 = 7?[4«^ — Kg + D(jk 0 — K 2 ) + «2 — K mArc] 

1 

71 D - 2n m 
Ak = K\ — K2 

D = 2jk 0 + 5 «: m - — . 

In the above equation, the subscripts n and t denote surface normal and tangential 
components of the field, U , t 2 are the principal surface directions and K\,k 2 are the 
corresponding curvatures, K m = (kj -\-k 2 )/2 is the mean curvature and K g = K\K 2 is the 
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Gaussian curvature. For /», = k, (spherical termination boundary). (17) reduces to the 
ABC derived by Webb and Kanellopoulos [57]. The latter authors have recent lv proposed 

a correction [58] for the implementation of these ABCs which should be incorporated in 
existing codes. 

Ipon substituting (17) into the surface integral of (G) and making use of vector 
differential and integral identities [59], we have 



• (H x n) dS = 


+ 


JJ s + a 2 (E^ ai f}dS 

J I x E scat ) 2 n dS 

JJs (V.E? Mf )[V-(7.E scat ) t )dS. 


(18) 


Next, expansion (6) or (14) is introduced and (18) is differentiated as in (7) to obtain 
the finite element equations (8) . Since the ABC is a local condition, the sparsity of 
the finite element matrix is preserved and the resulting system is again given by (12) 
with [Q] removed. However, in the case of arbitrarily curved boundaries, the matrix 
wi not be symmetric except for planar and spherical boundaries. The system will be 
also symmetric for cylindrical boundaries only if linear edge-based expansion functions 
are used. These ABCs have been extensively validated for a number of antenna and 
scattering configurations^], [60], Higher order ABCs have been proposed for a more 
effective suppression of the outgoing waves. These are often difficult to implement but 
as demonstrated by Senior et. al [61], they provide greater accuracy and effort being 
placed so much closer to the scatterer or radiator. ° 


3.3 Finite Element-Isotropic Absorber Method 

In accordance with this method, the outgoing waves are suppressed by an absorbing 
dielectric layer placed at some distance from the antenna as shown in Figure 6. Typ- 
ically, the layer is backed by metal and the finite element region extends all the way 
to the metal at which point the mesh is terminated by setting the tangential com- 
ponent of the total electric field to zero. This is equivalent to removing the integral 
over S 0 in (2)-(6). Obviously, the accuracy of the technique depends on how well the 
metal-backed absorber suppresses the incoming waves without introducing backward re- 
flections. The plane wave reflection coefficient of the absorber shown in Figure 6 has 
been minimized over the entire visible angular range [46], The fact that the dielectric 
has the same relative permittivity and permeability ensures that there is no reflection 
from the air-dielectric interface at normal incidence (perfect impedance match), but the 
performance degrades away from normal with the reflection coefficient reaching unity 
at grazing. This limitation led researchers to consider perfectly matched interfaces as 
discussed next. Nevertheless, even though the isotropic absorbers are not perfect, they 
are still useful in modeling antennas on doubly curved platforms. Figure 7(a) displays a 
sectoral microstrip patch printed on a conical surface and Figure 7(b) shows the antenna 
resonance behavior for different cone angles. We clearly observe that the resonance fre- 
quency drops with the cone angle for the same patch dimensions. It should be also 
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Figure 6: Illustration of an antenna terminated by an isotropic absorbing material 
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Figure 7: Sectoral microstrip patch 
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remarked that the computed resonance frequency is within 3.2 percent of that predicted 
h\' the approximate cavity model, and this is reasonable. 

In generating the results given in Figure 7. the computational domain was discretized 
using 2.358 prisms resulting in 3.790 degrees of freedom. One frequency run took 5.5 
minutes on an HP9000 (Model 715/04) workstation rated at 22 MFlops. 

3.4 Finite Element-Anisotropic Absorber Method 

This mesh termination method is similar to that described in the previous section except 
that the layer is comprised of an anisotropic lossy dielectric which has zero reflection 
coefficient at the air-dielectric interface over all incident angles and can therefore be 

considered as a perfectly matched interface (PMI) [20]. The intrinsic parameters of the 
PMI layer are 


a 


= tr 


a -j/3 0 

0 a- jfi 

0 0 


0 ' 

0 

i 

a-jP . 


where a (Z/Zq) 2 , Z is the intrinsic impedance of the medium being terminated by 
the absorber, whereas a, fi and the thickness of the layer are parameters which can be 
optimized for maximum absorption. 

The advantage of the anisotropic (over the isotropic) artificial absorber in terminating 
finite element meshes is illustrated in Figure 8 [62] where it is shown that the anisotropic 
absorber retains its low reflectivity at oblique incidences except near grazing (<f> = 0°). 
The performance shown in Figure 8 is based on a purely theoretical analysis but as 
shown in Figure 9, this performance can be realized in numerical simulations with careful 
choices of fit and the sampling rate N within the absorber. It has been found [63] that 
the reflectivity curve shown in Figure 9 is typical to most situations and the location 
of the minimum reflection coefficient can be predicted a priori An extensive numerical 
study based on a two dimensional planar absorber model demonstrated that fitj\ x (A r = 

A 0 / cos 4>) is a scalable parameter and that given the desired reflection coefficient Iffl 
the formulas [62] ’ 


Y = -0.0106|tf| + 0.0433 (19) 

N = 0.147e:rp[7.353y^] (20) 

can be used to choose fit/\ x and the minimum number of discrete samples to achieve 

the desired absorption. In these, \R\ is specified in dB and A c denotes the free space 

wavelength. These formulas are in agreement with the actual numerical results in Fig- 
ure 9 but it has yet to be determined how well they apply for curved perfectly matched 
layers which are placed conformal to scattering and radiating surfaces. Improvements 
to their absorptivity though can be attained by considering tapered layers and formulas 
similar to (19)-(20) are given by Legault [62] for one such tapered anisotropic absorber. 
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Anisotropic £ r = (T r = bT- (b - !/b) X x sin ^ 

(a) (b) 

Figure 8: Reflectivity of a A/4 thick planar metal-backed perfectly matched absorber 

(0 ~ ft) as a function of incidence angle, (a) Geometry, (b) plane wave reflection 
coefficient vs. angle 
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Figure 9. Reflection coefficient of the PM1 for terminating a microstrip line as extracted 
from a numerical implementation 
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4 Feed Modeling 


For scattering problems where the plane wave incidence is usuallv the ’source', the right - 
hand-side excitation has been explicitly given in (12) and will not be discussed further. 
However, for antenna parameter computations, the explicit form of {A' £ } in (11) will 
depend on the type of feeding scheme being employed. Below we discuss specific forms 
of {A e } corresponding to different feeding choices. 


Snnph Probe Feed: For thin substrates the coaxial cable feed may be simplified as a thin 

current filament of length / carrying an electric current II. Since this filament is located 
inside the cavity, the first term of the integral in (2) or (3) needs to be considered for 
this model. Specifically, the ith (global) entry of the excitation vector A', becomes 


A, - jk 0 Z 0 I 1 ■ W,(r), i = j u j 2t (21) 

where r is the location of the filament, m is the number of (non-metallic) element edges 
and j m is the global edge numbering index. In general, m such entries are associated 
with m element edges, and thus the index i goes from j x up to j m . This expression can 

be further reduced to A, = jk 0 Z 0 I l. provided that the ith edge is coincident with the 
current filament. 

Voltage Gap FeeeE This excitation is also referred to as a gap generator and amounts to 
specifying a priori^ the electric voltage V across the opening of the coax cable or any 
other gap. Since V = E • d, where d is a vector wdiose magnitude is the gap width, and 

E the electric field across the gap, we have that E, = — where cos6 t is equal to 1 

- dcosVi ^ 

l the i th edge is parallel to d. Numerically, this gap voltage model can be realized by 

first setting the diagonal term A,, of the [A] matrix equal to unity and the off-diagonal 
terms A tJ (i ^ j) to zero. For the right-hand-side vector, only the entry corresponding 
to the ith (global) edge across the gap is specified and set equal to the value E t whereas 
all other entries associated with edges not in the gap are set to zero. 

Coaxial Cable Fecd Model The simple probe feed of the coaxial cable is accurate only if 
the substrate is very thin. For thicker substrates, an improved feed model is necessary 
and this can be achieved by evaluating the functional 


F c = —jh 


° z °JL, 


(E xH )-zdS 


( 22 ) 


with 


E = — f. 


ho — 


see Figure 10) [64] 


H = —4>, 

T 

(23) 

yj^rc h 

rj e 0 H • 

Z 0 7 r 

(24) 
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F'gure 10: (a) Side view of a cavity-backed antenna with a coax cable feed; (b) Illustra- 
tion of the FEM mesh at the cavity-cable junction (the field is set to zero at the center 
conductor surface). 


In these expressions, c TC is the relative permittivity inside the cable, E and H are 
the electric and magnetic fields, respectively, measured at z = 0 and I 0 is the center 
conductor current. Also, (r,(f>,z) are the polar coordinates of a point in the cable with 
the center at r = 0. We observe that (24) is the desired constraint at the cable junction 
in terms of the new quantities h 0 and e 0 which can be used as new unknowns in place 
of the fields E and H. 

However, before introducing Fq into the system, it is necessary to relate e 0 and h 0 to 
the constant edge fields associated with the elements in the cavity region which border 
the cable aperture. Since the actual field has a 1 jr behavior in the cable, we find that 

AV = ~ a ) = e o^, * = N p (p =1,2, N c ) (25) 

where AV denotes the potential difference between the inner and outer surface of the 
cable and N p denotes the global number for the edge across the coax cable. When 
this condition is used in the functional Fc, it introduces the excitation into the finite 
element system without a need to extend the mesh inside the cable or to employ a 

fictitious current probe. The derivation of and its incorporation into the system is 

then a straightforward task [64]. As can be expected, the above feed model assumes the 
presence of only the dominant(TEM) mode at the cavity-cable junction, an assumption 
which may not be suitable for certain applications. Of course, the model can be improved 
by extending the mesh (say, a distance d) into the cable. The equi-potential condition 
will then be applied at z=-d, where all higher order modes vanish. 

Aperfure-Coup/ec? Microstrip Line Model: As shown in Figure 11, when the microstrip 
antenna is fed with a microstrip line network underneath the ground plane (cavity’s base) 
via a coupling aperture, special treatment of the feed structure must be considered in 
the FEM formulation. This is because the microstrip line is usually designed to have 
different size and shape as compared to the cavity’s geometries. Hence, the conventional 
simulation of treating the entire 3-D domain using a single type of element is not efficient 
or appropriate for this feed. 


20 


80 


Refining to Figure 11. it is appropriate to separate the computational domains be- 
cause of the small element size required in modeling the guided feed structure. As an 
example, let us consider a rectangular aperture which has been extensivelv employed in 
practice. The cavity fields may be discretized using tetrahedral elements, whereas in the 
microstnp line region rectangular bricks are the best candidate since the feed structure 
is rectangular in shape and the substrate is of constant thickness. Although both types 

Antenna fclments 


(0 

T 


St (II) 

i - — : 

Truncation Plane 

Coupling Aperture Sa 

Figure 11: Cross-section of an aperture coupled patch antenna, showing the cavity 
region I and the microstrip line region II for two different FEM computation domains. 

of elements employ edge-based field expansions, the meshes across the common area 
(coupling aperture) are different, and consequently some connectivity matrix must be 
introduced to relate the mesh edges across the aperture. However, since the aperture is 
very narrow, a ‘static’ field distribution may be assumed at any given frequency. There- 
fore, the potential concept may be again applied to relate the fields below and above 
the aperture. Specifically, the ‘equi-potentiaT continuity condition is enforced, and to 
demonstrate its implementation, let us first classify the slot edges as follows 

Tetrahedral Mesh (Cavity Mesh): 

• J = 1,2,3,... vertical edges 

• Ej 2 j — 1,2,3,... diagonal edges 

Brick Mesh (Feed Mesh): 

• Ej j = 1,2,3,... vertical edges only 
Then the ‘equi-potential’ continuity condition requires that 

E f = tjE' 

E T = Yd^ + t ^EU) 

in which 


21 


81 


whereas f and d are the lengths of the vertical and diagonal edges, respectively. That 
is, t is simply the width of the narrow rectangular aperture. The coefficient ij is equal 
to ±1 depending on the sign conventions associated with the meshes above and below 
the coupling aperture. 

This connects ity scheme can be extended to entirely different computational do- 
mains. Also, it is apparent that the approach makes the FEM implementation straight- 
forward for different geometry/size domains that would be significantly inefficient if only 
one type of elements were used for modeling the structure. In addition, the technique 
ensures a good system condition since the number of distorted elements in the mesh is 
minimized. 


5 Parallelization 

When considering 3D problems of practical interest, the unknown count of the compu- 
tational domain can easily reach several million degrees of freedom. The sparsity of the 
FEM system (particularly for the FE-ABC and FE-AA methods) makes possible the 
storage of such large scale problems but even at 0(N) computational demands, their 
practical solution requires efficient use of parallel and vector platforms. Modern com- 
puting platforms can now deliver sustained speeds of several GFlops and CPU speeds 
m the T flops range are within our reach. The inherent sparse matrices of PDE meth- 
ods are particularly suited for execution on multiprocessor and vector platforms but 
the exploitation of these processors requires special storage schemes and techniques to 
perform the matrix-vector product required in the iterative algorithms at the Flop rates 
sustained on these multiprocessors. 

To parallelize and vectorize the FEM codes, it is essential to first optimize the ex- 
ecution of the iterative solvers which typically take-up 90% of the CPU time. Among 
them, the conjugate gradient algorithms (CG, BCG, CGS and QMR) have been found 
very attractive and a brief comparison of the pros and cons for these is given in [45]. The 
Generalized Minimal Residual Method (GMRES) is another iterative solver which can 
exhibit faster convergence rates. However, it stores the direction vectors and as a result 
it requires much higher storage. For the discussion below we will primarily concentrate 
on the BCG and QMR algorithms and we note that the symmetric form of BCG re- 
quires minimal number of arithmetic operations (see Table 2). A disadvantage of the 
BCG is its erratic convergence pattern whereas the QMR has smooth and monotonic 
convergence. However, neither BCG nor QMR can guarantee convergence and typically 
they both converge or not for the same problem. When considering the parallelization 
of a completely sparse system such as that resulting from the FE-ABC method, the 
following issues must be addressed: 

Storage^ of Sparse System: The performance of the code is strongly dependent on the 
employed storage scheme. Since a typical FEM matrix has about 8.5 N ee or so non-zero 
entries, it is essential that the non-zero elements be stored in a manner that keeps the 
storage requirements nearly equal to the non-zero entries and minimizes inter-processor 
communications. The ITPACK [65] and the compressed row storage (CRS) schemes are 
most appropriate for parallel computing. The ITPACK storage format is most efficient 
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Complex 

Operation 

★ 

+ 

Matrix-vector Products 

nze 

nzc-X 

j Vector Updates 

IN 

3X 

Dot Products 

3N 

3X 

Total # of Operations 

29 X 

3X 


X = # of unknowns 

nze = # of nonzero matrix elements 

Table 2: Floating Point Operations of BCG Per Iteration 


for retrieving the matrix elements and certainly the choice method when the number 
of non- zero elements are nearly equal for every matrix row. Basically, the ITPACK 
format casts the FEM matrix in a smaller rectangular matrix having the same rows as 
the original matrix and can be unpacked by referring to a pointer integer matrix of the 
same size. However, this rectangular matrix can contain as much as 50% zeros which 
results in space wastage. By using a modified ITPACK scheme, space wastage can be 
reduced down to 30%. Even with less wastage, the CRS format may be the most efficient 
storage scheme with some compromise in CPU speed. It amounts to storing [.4] as a 
single long row which can be uncompressed using two integer pointer arrays. For the 
symmetric BCG algorithm, the CRS format results in only 8.5 N complex numbers and 
N integers. However, it should be pointed out that the CRS format is not appropriate for 
vector processors such as the C-90. For vectorization, it is best to organize the storage in 
sections of long vectors and to achieve this for our type of matrices the jagged diagonal 
format [66] appears to work best. Using this format the rows are reordered so that the 
row with the maximum number of non- zeros is placed at the top of the matrix and 
rows with the least non-zero entries are shuffled to the bottom. This reordering greatly 
enhances vectorization because it allows tiling of the shorter rows to yield very long 
vector lengths in the matrix-vector multiplication phase. Specifically, for some problem 
the jagged diagonal storage format allowed the matrix-vector multiplication routine to 
run at about 275 MFlops on a Cray C-90 whereas the same routine ran at 60 MFlops 
using the CRS format. The dot product speeds and the vector updates reached 550 
MFlops and 600 MFlops for the same problem. Table 3 provides a relative comparison 
of CPU estimates on various computers. 

/niter-processor CommtimcaU'ons: For distributed memory platforms, the method of par- 
titioning the stiffness matrix [A] among the processors, the chosen storage scheme and 
the inherent unstructured sparsity of [.4] are all crucial to the overall speed of the code. 
An approach that has worked well on massively parallel processors (such as the SP-2, 
Intel Paragon, Convex Exemplar) is that of assigning each processor a section of the 
matrix and by dividing the vectors among the P processors. Thus, each processor is 
responsible for carrying out the matrix-vector product for the block of the matrix it 
owns. However, the iterate vector is subdivided among all processors, and therefore 
narrow-band or structured sparse matrices have an advantage because they reduce in- 
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Processors 

# of Processors. P 

Tieni f //-secs/it oral ion/unknown ) f 

Cray C-90 

1 (275 MFlops) 

0.55 

KSR 1 

28 

1.28 

58 

0.57 

Intel Paragon 

8 ^ 

3.42 

16 

1.99 

32 

1.38 

IBM SP-1 

4 

1.47 


Table 3: CPU Time Per Unknown for Solving Typical FE-ABC Systems 


terprocessor communication. Since typical FEM matrices are unstructured, algorithms 
such as the Recursive Spectral Bisection (RSB) have been found very effective in re- 
ducing inter-processor communication. However, the standard Gibbs-Pool-Stockmeyer 
profile reduction algorithm has been found even more effective in reducing the initial 
FE-ABC matrix (see Figure 3) to banded form as illustrated in Figure 12. This type of 
matrix reordering can deliver speed-ups as close to linear as possible. 



Figure 12: Reduced bandwidth of the FE-ABC system after application of the Gibbs- 
Pool-Stockmeyer profile reduction algorithm 

Matrix Preconditioning: Preconditioned iterative solvers are intended to improve the 
convergence rate of the algorithm. At times, preconditioners are necessary as may 
be the case with some dielectrically loaded structures. However, for relatively small 
systems (less than 100,000 unknowns) it has been found that diagonal preconditioning 
is typically most effective and should always be applied. This preconditioning amounts 
to normalizing each row by the largest element, but even this simple operation can lead to 
substantial convergence speed-ups. Block and incomplete LU preconditioners are more 
effective in improving the convergence of the solver but are more costly to implement 
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and oih* must judge on the overall CPU requirements rather than on the improved 
convergence alone. For example, the incomplete LU preconditioner given in [67] reduced 
the iterations to 1/3 of those needed with diagonal preconditioning. However, each 
iteration was 3 times more expensive due to the triangular solver bottleneck. 

6 Additional Applications 

We choose three more examples to demonstrate the capability of the hybrid finite element 
methods. 

Scattering by a Large Cone-Sphere: A cone-sphere is basically a hemisphere at- 
tached to a cone. This is a difficult geometry to mesh since a surface singularity exists at 
the tip of the cone. The singularity can be removed in two ways: i) by creating a small 
region near the tip and detaching it from the surface or ii) by chopping off a small part 
near the tip of the cone. The second option inevitably leads to small inaccuracies for 
backscatter from the conical tip; however, we chose this option since the conical angle 
in our tested geometry was extremely small (around 7°) and the mesh generator failed 
to mesh the first case on numerous occasions. In Figure 13, we plot the backscatter pat- 
terns of a 4. 5 A long cone-sphere having a radius of 0.5A for 66 polarization. The mesh 
truncation surface is a rectangular box placed 0.4A from the surface of the cone-sphere. 
As seen, the far-field results compare extremely well with computations from a body of 
revolution code [68]. 

Frequency Selective Surfaces (FSS): FSS structures [69] are arrays of tightly packed 
periodic elements which are typically sandwiched between dielectric layers. The periodic 
elements may be of printed form or slot configurations designed to resonate at specific 
frequencies. As such, they are penetrable around the element resonances and become 
completely reflecting at other frequencies. To meet bandwidth design specifications, 
stacked element arrays may be used in conjunction with dielectric layer loading. Here 
we consider the analysis of FSS structures via the FE-BI method. Because of the fine 
geometrical detail associated with the FSS surface, the finite element method has yet 
to be applied for the characterization of FSS structures, but use of prismatic elements 
makes this a much easier task. Of particular interest in FSS design is the determination 
of the transmission coefficient as a function of frequency, and since the array is periodic, 
it suffices to consider a single cell of the FSS. For computing the transmission coefficient 
T i the periodic cell is placed in a cavity as shown in Figure 14 and the structure is excited 
by a plane wave impinging at normal incidence. Assuming that near resonance the wave 
transmitted through the FSS screen will retain its TEM character, the transmission 
coefficient of the FSS panel can be approximated as 

t2$ = 10 log — 
a 

where a is the reflection coefficient of the absorber placed at the bottom of the cavity 
and should be kept small (< 0.1) to suppress higher order interactions. By adding the 
next higher order interaction, a more accurate expression for the transmission coefficient 
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Figure 13: Backscatter pattern of a perfectly conducting conesphere for <f><f> and 00 
polarizations. Black dots indicate computed values using the FE-ABC code (referred 
to as FEMATS) and the solid line represents data from a body of revolution code [68]. 
Mesh termination surface is a rectangular box. 
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Figure 14: Upper figure: geometry of the multilayer frequency selective surface (FSS) 
used for modeling; lower figure: measured and calculated transmission coefficient 

through the FSS structure 


is 


TdB ~ Tjg + 10 log [l - cr(l - T^)] . 

The above FSS modeling approach was applied for a characterization of multi-layered 
slot FSS structures. The geometry of the multilayer radome is given in Fig. 14. The 
total thickness of the FSS was 6.3072cm and is comprised of two slot arrays (of the same 
geometry) sandwiched within the dielectric layers. For modeling purpose, a 1.54cm thick 
absorber is placed below the FSS as shown in Fig 14. It is seen that the results generated 
by the FE-BI method are in good agreement with the measurements [70]. 

Scattering by Jet Engine Inlets: 

This example demonstrates a hybridization of the finite element and high frequency 
methods for scattering by jet engine structures[71]. In this hybridization, the jet en- 
gine itself is modeled by the FEM taking advantage of the geometrical adaptability of 



finite methods. The fields generated by the FEM are then propagated to the mouth 
of the jet engine using a raj' technique such as the geometrical theory of diffraction or 
the shooting and bouncing ray methods. To analyze the scattering by the jet engine, 
the blade structure is placed in a computational domain which is truncated by a per- 
fectly matched absorbers at about 0.5 wavelengths from the engine face. By exciting 
the engine on a mode by mode basis, the appropriate modal scattering matrix is gen- 
erated using the procedure described in [71]. However, in this case the computational 
domain is confined to a single blade lobe (cylindrical periodic cell) using appropriate 
phase boundary conditions at the boundaries of the periodic cells. This reduces the 
computational requirements to very low orders and permits the simulation of practi- 
cal jet engine geometries. Given the modal scattering matrix of the jet engine, high 
frequency techniques can be readily used for propagating the fields to the inlet open- 
ing as described by Karty et.al.[72]. A demonstration on the accuracy of this hybrid 
technique is given in Figure 15 where we show the computed scattered field due to a 
single mode excitation. The specific engine-blade geometry was placed in a cylindrical 
metal-backed cavity as shown in Figure 15. It consisted of 20 straight blades mounted 
radiallly on a cylindrical hub 9.25” in radius with the blades extending all the way to 
the outer wall of the cavity which was 14.375” in radius. The blades were 2 degrees 
thick and 2.75 v long, whereastheentireinletdepthwas\9.685 n and we remark that this 
configuration is similar to an F4 fighter engine. The shown computation was carried 
out at 4.348GHz and involved the solution of an FEM system with only 87,000 edges 
(64,000 elements) since a single periodic lobe was descretized but the overall structure 
spans about 4.5 wavelengths in diameter. As seen the agreement between the FEM com- 
putations and those of the reference modal solution are in excellent agreement, nearly 
overlaying each other. 

7 Conclusion 

In this paper we reviewed hybrid finite element methods as applied to electromagnetic 
scattering and radiation problems. Much of the emphasis dealt with the various mesh 
truncations schemes and we presented an up-to-date account of these schemes. The 
usual finite element-boundary integral method was presented and new developments for 
reducing the CPU requirements of this technique using the fast integral methods were 
discussed. Antenna feed modeling in the context of the finite element method had not 
been discussed before and for the first time we presented an overview of the model- 
ing approaches for the most popular antenna feeds, including aperture coupled feeds. 
Parallelization will continue to play an increasingly greater role and a section was in- 
cluded discussing our experiences for better implementation of finite element codes on 
distributed and vector architectures. A number of examples illustrating the success- 
ful application of the finite element method were included throughout the paper and 
these were intended to demonstrate the method’s geometrical adaptability and inherent 
capability to treat highly heterogeneous structures. 

As can be expected, issues relating to mesh truncation, mixing of elements [73], do- 
main decomposition[74, 75], robustness, adaptive refinement [76], accuracy, error control, 
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Figure 15: Scattering by a cylindrical guide terminated with a straight bladed engine 
configuration. 
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eed modeling and parallelization for large scale simulations will continue to dominate 
future research and developments relating to partial differential ec|uation methods. An 
apparent advantage of the finite element method is its potential hvhridization with all 
other frequency domain methods. Future applications of the finite element method are 
) e } to make greater use of hybridization techniques aimed at increasing the method's 

accuracy and efficiency while retaining its inherent geometrical adaptabilitv and ease in 
handling materials. 
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